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Abstract 

An exact analysis of the N-dimensional blackbody radiation process in a 
nonextensive a la Tsallis scenario is performed for values of the nonextensive's 
index in the range (0 < q < 1). The recently advanced "Optimal Lagrange 
Multipliers" (OLM) technique has been employed. The results are consistent 
with those of the extensive, q = 1 case. The generalization of the celebrated 
laws of Planck, Stefan-Boltzmann, and Wien are investigated. 
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I. INTRODUCTION 



Planck's blackbody radiation studies constitute one of the milestones in the history of 
Physics. Planck's law is satisfactorily accounted by recourse of Bose-Einstein statistics. The 
possible existence of small deviations from this law in the cosmic blackbody radiation has 
recently been discussed. They would have arisen, at the time of the matter-radiation de- 
coupling |IJ as a consequence of the fact that long-range interactions seem to be associated to 
a nonextensive scenario PH. The concomitant thermostatistical treatment [Q-[nj is by now 
recognized as a new paradigm for statistical mechanical considerations. It revolves around 
the concept of Tsallis' information measure S, |, a generalization of Shannon's one, that 
depends upon a real index q and becomes identical to Shannon's measure for the particular 
value q = 1. 

The blackbody radiation problem was discussed within the Tsallis' nonextensive frame- 
work in [0. It was there conjectured that the deviations from Planck's law detected in the 
pertinent data could be attributed to a long-range gravitational influence. The study 
reported in |J employs the so-called Curado-Tsallis unnormalized expectation values M in 



the limit q — > 1. Afterwards, Lenzi |12| advanced an exact treatment of the problem, also 
within the Curado-Tsallis framework. 

Nowadays, it is believed that this framework has been superseded by the so-called nor- 
malized one, advanced in [13,14 , which seems to exhibit important advantages ||. This 
normalized treatment, in turn, has been considerably improved by the so-called "Optimal 
Lagrange Multipliers" (OLM) approach ||15|| . It is then natural to revisit the problem from 
such a new context. 

In such a vein we shall consider blackbody radiation in equilibrium within an enclosure of 
volume V with the goal of ascertaining the possible g-dependence of the Planck spectrum, 
Steffan-Boltzmann's law and Wien's one. The OLM treatment seems to be indicated in 
view of the findings of fl6lJT7|| regarding the particular nature of the Lagrange multiplier 
associated to the temperature. 
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The paper is organized as follows: In Sect. [TJ we present a brief OLM primer, while in 
Sect. Ill we obtain the partition function, internal energy and energy density of our problem. 
In Sect. IV first-order corrections in 1 — q are discussed and, finally, some conclusions are 
drawn in Sect. |V[ 



II. MAIN RESULTS OF THE OLM FORMALISM 



For a most general quantal treatment, in a basis-independent way, consideration is re- 
quired of the density operator p that maximizes Tsallis' entropy 

S q 1 - Tr(p«) 



(1) 

M) denote the 



k q — 1 

subject to the M generalized expectation values {Oj)-, where Oj (j = 1, . 
relevant observables. 

Tsallis' normalized probability distribution |H| , is obtained by following the well known 



MaxEnt route [18]. Instead of effecting the variational treatment of Tsallis- Mendes-Plastino 



(TMP) ||13|| , we pursue the alternative path developed in [15|. One maximizes Tsallis' 
generalized entropy (P PJ5|JT9|] subject to the constraints ||{n 



Tr{p) = 1 



Tr 



whose generalized expectation values 



O, 



p q [o 3 -(o. 



0. 



(2) 
(3) 



Tr(p q Oj 



<? Tr(p q ) 

are (assumedly) a priori known. The resulting density operator reads 



(4) 



p = K 1 



M , N 

l-(l-?)£^<5;-(6;)J 



l-q 



(5) 



where {Xj} stands for the Optimal Lagrange Multipliers' set and Z q is the partition function 



Z q = Tr 



M 



i-(i-ff)E^ i(Oi-(o. 



1-9 



(6) 



3 



It is shown in ITL5H that 



Tr(p q 



and that Tsallis' entropy can be cast as 



(7) 



S q = k \n q Z q , 



(8) 



with ln q Z q = (1 — Zz q )/{q — 1). These results coincide with those of TMP |13| in their 



normalized treatment. If, following |Tj|, we define now 



\n q Z q = \n q Z q -Zl- q Y,^ (Oj 



(9) 



we are straightforwardly led to [15 







S r , 



did,) \k 



(In, Z q 



d 
d\ 3 



■z±- q (d 



(10) 

(ii) 



Equations (P~0| ) and ( O ) are modified Information Theory relations like those used to 
build up an a la Jaynes Statistical Mechanics. The basic Legendre-structure relations 
can be recovered in the limit q — > 1. 



III. BLACKBODY RADIATION 



A. The standard situation q = 1 



In order to analyze the blackbody radiation in the generalized statistical context we 
first consider the standard situation, q — 1. We look for the equilibrium properties of the 
blackbody electromagnetic radiation. The appropriate thermodynamical variables are the 



volume V and the temperature T [pO| , pT 



The Hamiltonian of the electromagnetic field, in which there are rik, e photons of momen- 
tum k and polarization e, is given by 



W = X)7iam fe>e , (12) 

k,e 

where the frequency is lj = c |k| (c is light's speed) and n^^ = 0, 1, 2, . . ., with no restrictions 
on {n^} . The partition function reads 

Z 1 = Tr(e^). (13) 



For a macroscopic volume we define a density of states in the n-dimensional space g n (uj) = 
A n u) n - X with 

A 2TnV (U) 

n (47rc 2 )"/ 2 r(n/2)' { ' 

where r n = n — 1 is the number of linear-independent polarizations. Eq. fll3|) can be written 



as 



Z 1 = exp|y dujg n {u) In [1 — exp (— /3?kj)]j = e^ n , (15) 



where 



T A 



with 



B. Generalized situation q ^ 1 (0 < q < 1) 



1. Generalized Partition Function 



The OLM-Tsallis generalized partition function (see Section |T|) is 



Z q = Tr 



(16) 



I n = - J™ dxx n ~ l In (l - e- x ) = r(n)C(n + 1), (17) 
C stands for the Riemann z-function and V is the Gamma function. 



l-(l- q )p(H- Uq )y-\ (18) 



where H is the Hamiltonian fll2]) , U q is the mean energy (§]), and /3 is the optimal multiplier 
Lagrange. The photons' chemical potential is, of course, taken as zero. 



With the aim of calculating Z q we follow the steps of [12]. We work again in the large 



volume limit and regard Eq. (|T8|) as an integral that can be evaluated using the relation [22 



e ab b i-z °? itb i for 6 > 

d h — z = i ( 19 ) 

for6<0, 



2n J (a + it)" 



with a > 0, Re(z) > 0, and -n/2 < arg(a + it) < ir/2. 

If we set o = l, 6=1 — (1 — q)(3{H — U q ) (the cut-off condition is naturally fulfilled [Q) 
and z — 1/(1 — q) + 1, the generalized partition function adopts the appearance 

oo 

Z,(/3) = | dtK q (t)e^ u "Z 1 0), (20) 

— oo 

with Z x given by Eq. flT5|), 

^ m _ r[(2-g)/(i-g)]ex P (i + ^) r9n 
elj 2tt(i + ^)( 2 -'?)/( 1 ^) ' 1 j 

and 

^=(l+it)(l-g)/3. (22) 
In order to evaluate the integral in ( p0|) we expand the exponential term and obtain 



= £ g rip -,)/(!-,)] J ^'■+"-~» , (23) 

— 77? I /7T / i'i i • « \ ^ — -+nm 



where £ n is given by Eq. Qlq) . 

Using again (|l^) with 6 = 1 + (1 — q)/3U q , a = 1, and z = (2 — g)/(l — g) + ?im we arrive 

to 

= E ^ rrr9 r[(2 :' )/(1 7' )] 1 [1 + (1 - g)/?^]^ +nm , (24) 
n^o ™ r K 2 - <?)/(! - 9) + nm\ 

and notice that an additional cut-off's like condition must be considered 1 + (1 — q)/3U q > 0. 

Otherwise, Z q (f3) = 0. 



2. The internal energy 



The generalized internal energy can be cast in the fashion 



<? rp 1 



(25) 



where we have introduced the definitions 



T\=Tr \ l-(l- q )p(H-U q ) 



Tr 



l-(l-q)P{H-U q 
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H 



(26) 
(27) 



We first analyze T 2 . The expression 1 — (1 — q)f3{H — U q ) 1 q is evaluated by recourse 
to fll9|) , where we set a = 1, 6 = 1 — (1 — q)f3(H — U q ), and z = 1/(1 — g). We are thus led to 

^_r[l/(l-g)] 



l-(l- g )/3(#-^ 
The trace operation yields afterwards 

r[i/(i-g)] 



2tt 



X e (l+tt)[l-(l-q)f3(H-U q )] 

dt i 



;i + a) 



(28) 
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2vr 



e (l+ii)[l+(l-g)/3C/ 9 ] 



^l(/3), 



(29) 



{1+it) — 

with /3 defined by Eq. Q2"^). Expanding exponents and employing once again (|IP|) leads 
finally to 



m=0 ml T[l/(1 -g) +nm 
As for T 2 , we recast Tj as 



(30) 



r[l/(l-g)] 7 e (i+it)[i+(i-<i)PU q ] 



27T _^ (i + it )~ q 
We proceed as follows: i) take advantage of the fact that 



Tr 



Z 1 0)H 



(31) 



Tr 



Zi(j3)H 



d(3 



(32) 



ii) use Eq. (|T5|) , and iii) expand again the exponential term. We find 



_ nr[l/(l-g)] ^c +1 
1 " 2tt(1 - g)/3 ^ m! _ 



dt- 



3 (l+ii)[l+(l-?)/3E/g] 



[1 + it) 



(33) 



The integral is evaluated as before, which yields 

Tl = (T^ |„^ r [V(1 -T^i) + i] 11 + (1 " ■ « 34 » 

so that the internal energy becomes 

V SI 1 [1 + (\ _ ,/|tf/7 in(m+l)+l 

n m\ r[l/(l-g)+n(m+l)+l] I 1 ^ \ L H)t Ju q\ 
jj _ '"an m=0 L / v h;-t \ (35) 

m=0 

although compliance with the Tsallis'cut-off condition 1 + (1 — q)0U q > is always to be 
demanded. The above is a non-linear equation to be tackled in numerical fashion. Some 
pertinent results are displayed in Fig. [I], for different g-values, they include the standard 
q — 1 case that yields the Steffan Boltzmann law 

U x = ctT 4 (36) 

with a = 5.67 x 10~ 8 W/m 2 K A , the so-called Steffan- Boltzmann constant. 

We notice that a Steffan Boltzmann-like law holds, if kT is small enough, for all q G 
(0 < q < 1), with a = o~(q). For high T- values and q 1 the power law behavior seems to 
persist, but the corresponding power, let us call it a q , is no longer equal to 4. Fig. || is a 
graph of a q versus q for 

U q oc T a ". (37) 

Fig. ^] (a) shows a q for small fcT-values. The validity of Stefan Boltzmann's Law for 
a wide range of g-values is easily appreciated. Fig. |2| (b) depicts, instead, a q for large 
kT- values. In this limit, important deviations from the a q = 4 value are detected. 

An intermediate temperature range is observed in Fig. [I] in which the power law behavior 
(of U q vs. kT) is violated. The corresponding transition is the more abrupt the larger \q— 1|. 
As will be seen in the next section, these violations can be attributed to a Tsallis cut-off in 
the energy densities. 
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3. Energy densities 



The generalized spectral energy distribution of blackbody radiation u q is given by 

roc 

U q = / doou q . (38) 





In order to obtain u q we work with U q = T1/T2, with Ti and T2 given by Eqs. fl26|) and 
(p7|), respectively. The denominator will not suffer any change, but in the numerator we will 
keep the integral form of Z\ as given in Eq. (|l5l). Then, Eq. fl32| ) will lead to 

Tr (e^tf) = ^ / ^ B 7— ^. (39) 



allowing us to cast T\ in the fashion 

V\\l(\-a\\ 7 7 (i+it)[i-(i-q)/3(hu-u q )] £, n 
T = MW gJJ M A A ^ (4Q) 



Recourse to the identity 



1 00 

E (^^) , (41) 



1 - s=0 
yields now 

= MVU gJJ ^ /■ Vk eft- j . (42) 

If we use Eq. (|19]) the integral above is easily calculated, and we obtain 







According to (38) the ratio ^ implies that 



u q = hA n u n ^ , (44) 



where 



m=0 



1 Tl II! 



o C [1 + (1 - g)/^ 

m m! r[l/(l - ?) + nm] ' 1 J 



8=1 



(l-q)f3[hm(s + l)-U q ] 
1 + (1 - g)/?[/ 9 



t-2 — \-nm 

l-q 



(46) 



Figs. and [| are plots of u q vs. the frequency uj. Notice that they quite resemble the 
semi-empirical Planck law for small \q — l|-values. Fig. [5] clearly shows that the frequency 
shifts of the maxima of the above referred to curves (located at uj — uj*) comply with the 
Wien law 

uj* oc T, (47) 

for such g-values. Up to q — 0.8 this nice behaviour persists. As q departs from 1 an 
anomalous behaviour ensues. Local maxima begin to appear in the curves u q (uj) (Figs. |3| 
and |j). As T grows the global maximum is shifted and the secondary ones become more 
pronounced, until one of them replaces the former one as the maximum maximorum. An 
abrupt jump in the curve uj* vs. T ensues as a consequence (Fig. ||). 

The departures from Planck's law as \q — 1| grows, could in principle be subject to 
experimental verification, thus q would be fixed. 

The above referred to discontinuities in the "Wien-like" behaviour are due to the Tsallis' 
cut-off operative in terms of the sums T\ and T 2 . 



IV. FIRST ORDER CORRECTIONS 



We have seen in the preceding section that (|44| ) yields, for g-values close to unity, results 
that quite resemble the ones given by Planck's law. It is then reasonable to look for a 
perturbative expansion in 1 — q. 

Let us start with S m (Cf. Eq. ([|6|)) and write it in the form 

4- s\ 1 1 

(48) 



=0 



. f3huo{l + s) 



1 + (1 - q)(3U q 



A first order expansion of the logarithm yields 



s ^lr{- [q+nm{1 - q)] T^M' m 

a power series in s of guaranteed convergence for 1 — q < 1. Neglecting the term nm(l — q) 
on account of the fact that the concomitant series rapidly converges so that terms for which 
m could compete with 1/(1 — q) are indeed negligible, we have 
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and replacement in Eq. (pHD yields 



HA n u n 



u q m 



exp 



qfthu) 



- 1 



, (51) 

l+{\-q)pU q \ ' 

a first order correction for a "generalized Planck law" that yields the classical result for 
q — > 1. Eq. ( pi]) provides one then with an approximate energy density. 

Using now this equation as a starting point we also get first order corrections in 1 — q to 
the Wien law. We can locate the maximum with respect to uj of Eq. (O) with the auxiliary 
definition 

uqf3Ti 



x 



and immediately find 



l-(l-q)/3U q 



e - + - = l, 
n 



(52) 



(53) 



whose solution is a constant b for each fixed value of n. For n = 3 we find, for instance, 
b = 2.82. The maximum of u q (u) for different T's is located at distinct frequencies according 
to 



bk 

0Ji = —Ti 
qh 



l-qU q 



(54) 



q h 

a first order correction to Wien's law. 

If we integrate Eq. fl5"T|) over frequencies we find U q as a function of T, and a concomitant 
first order correction to Planck's law 

n+l 



U q = hA n a 



1 + (1 - q)0U q 
qph 



(55) 



that reduces trivially to Planck's one in the limit q — > 1. We call a the integral 

a= dx-— = T(n + l)((n + l), (56) 
Jo e x — 1 

with x given by (|52|). 



It may be convenient to point out that the first order relationships here obtained ignore 
possible cut-off problems. Anyway, in the < q < 1 range (with q — > 1) that is relevant to 
the present considerations, such problems should not arise. 

V. CONCLUSIONS 

We have presented here exact results for the mean energy and the energy density of 
the blackbody problem in a nonextensive environment. A Stefan Boltzmann-like law is 
obtained, for all < q < 1, although with a g-dependent SB constant cr(q), for almost the 
whole temperature range, although for some temperatures a deviation from the power law 
behaviour is detected. 

The energy density curves correspond to Planck's law for q values close to unity. For other 
g's a gradual departure from the typical Planck curve is observed, which can be attributed 
to Tsallis' cut-off effects. Wien's law is also obtained for q values close to unity. 

Finally, we have presented first order perturbative results in the limit q — > 1 for the 
energy density u q (the parameter is, of course, 1 — q). These could be of some utility in 
obtaining first order nonextensive corrections to the Bose-Einstein equation. 
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FIGURES 




kT 

FIG. 1. The internal energy Uq as a function of T = 1/(3 for different values of the 
non-extensivity parameter q. 




FIG. 2. Power law exponent a q as a function of q (see text for details). The graph (a) corre- 
sponds to small temperatures, and the (b)-one to high T-values. 
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to (0 (0 



FIG. 4. Energy density u q as a function of the frequency uj. Each plot correspond to different 
values of q parameter. In all of them, kT = 10 2 was kept fixed. 
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q=1.0 




kT 

FIG. 5. We depict the frequency oj* for which the energy density reaches its maximum, as a 
function of kT, for different values of q. 
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